Generation of mechanical squeezing via magnetic dipoles on cantilevers 
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A scheme to squeeze the center-of-mass motional quadratures of a quantum mechanical oscillator 
below its standard quantum limit is proposed and analyzed theoretically. It relies on the dipole- 
dipole coupling between a magnetic dipole mounted on the tip of a cantilever to equally oriented 
dipoles located on a mesoscopic tuning fork. We also investigate the influence of several sources 
of noise on the achievable squeezing, including classical noise in the driving fork and the clamping 
noise in the oscillator. A detection of the state of the cantilever based on state transfer to a light 
field is considered. We investigate possible limitations of that scheme. 



I. INTRODUCTION 



Mechanical cantilevers have a long and rich history as 
force and field meters. In recent decades, microf abdica- 
tion and nanotechnology have resulted in spectacular ad- 
vances with the development of tip microscopy [l[ , result- 
ing in the measurement of feeble forces and the imagingof 
single atoms, nanoscale magnetic resonance imaging |2|, 
and single spin detection |3[, to mention just a few ap- 
plications. With recent progress in cavity optomechanics 
and the successful cooling of at least three micromechan- 
ical systems to the deep quantum regime with just 
a fraction of a phonon of center-of-mass excitation left, 
mechanical sensing is at the threshold of an important 
new breakthrough. 

One potential advantage of operating micromechani- 
cal sensors in the quantum regime is that this opens the 
way to measurement techniques that can circumvent the 
standard quantum limit. These techniques rely on the 
capability to locate the unavoidable quantum noise in a 
quadrature of the field to be measured that does not in- 
teract with the measuring apparatus, and to implement 
back-action evading techniques that prevent that noise to 
feed back into the outcome of successive measurements. 
In quantum optics, the most famous states that permit to 
achieve this goal are single-mode squeezed states, where 
the variance of one of the quadratures X = (a + at)/ 2 
or X2 = (a — at)/2i of the light field is below 1/4, 
with (AAfAAf) > 1/16. Here a and a) are normal- 
ized bosonic annihilation and creation operators of the 
field mode, and (AX?) = (Xf) - (X,) 2 . 

While the classical noise squeezing of micromechani- 
cal oscillators has already been achieved @, [1| , squeezing 
below the standard quantum limit still remains to be 
demonstrated. A number of techniques have been pro- 
posed, including conditional squeezing using a paramet- 
rically coupled electromagnetic cavity driven by one Q 
or two sidebands [l(J E2 detuned from the cavity reso- 
nance by the mechanical oscillation frequency, and are 
expected to be demonstrated experimentally in the near 
future. This paper proposes an alternative scheme where 
a nanoscale cantilever can be prepared in a squeezed state 
by purely mechanical means via the nonlincarity provided 




FIG. 1. Classical tuning fork magnetically coupled to a 
quantum-mechanical cantilever. The equilibrium position of 
the cantilever is equidistant from the two extremities of the 
nanoscale tuning fork. 



by the magnetic dipole-dipole interaction between a vi- 
brating classical fork and the cantilever. We also show 
that the squeezing can be detected via state transfer to 
an optical field coupled to the cantilever in an optical res- 
onator configuration fl2l - [l4l ]. We discuss the impact of 
the various sources of noise, including both the noise of 
the classical driving force and the clamping noise of the 
cantilever and also comment on limitations to the state 
transfer scheme resulting from the dynamics of the light 
field and from the opening of an associated dissipation 
channel. 



The remainder of this paper is organized as follows: 
Section II introduces our model system and demonstrates 
how the anharmonic potential that describes the mag- 
netic coupling between the classical fork and the quantum 
mechanical cantilever results in quadrature squeezing of 
the cantilever motion under appropriate conditions. Sec- 
tion III analyzes the robustness of the system against 
various sources of noise. Section IV discusses the beam- 
splitter state transfer mechanism, and section V is a con- 
clusion and outlook. Calculational details are relegated 
to two appendices. 
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II. MODEL SYSTEM 

We consider a nanomechanical system consisting of a 
cantilever magnetically coupled to a classical nanoscalc 
tuning fork, as shown in Fig. [TJ the coupling being real- 
ized via point-like magnetic dipolcs located at both ex- 
tremities of the fork as well as on the cantilever. The 
oscillation direction of the cantilever and the fork are 
taken to be the x-axis and z-axis, respectively, and the 
magnetic dipoles are assumed to point in the positive 
z-direction. We consider specifically the center-of-mass 
mode of vibration of the fork, in which the distance be- 
tween its two extremities remains constant. This ar- 
rangement provides a stable mode of operation for ex- 
perimentally reasonable parameters. 

Denoting the displacement of the cantilever from its 
equilibrium position by x c , the interaction energy be- 
tween the three magnetic dipoles can be written as 
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where d / and d c are the dipoles moments of the magnets 
attached to the tips of the fork and to the cantilever, 
respectively, and £± are the distances between the two 
tips of the fork and the cantilever along the ^-direction. 



A. Symmetric case 

We assume first that the equilibrium position of the 
nanomagnet on the cantilever is equidistant from the two 
tips of the fork. For that setup the attractive forces from 
the two dipoles acting on the cantilever cancel each other, 
but for any departure from that situation this is no longer 
the case, so that the stability of the system must be en- 
forced by the stiffness of the cantilever. 

For the mode of vibration under consideration we have 
in the symmetric case 
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where Zf denotes the displacement of the tuning fork tips 
from their equilibrium position and £ = (£+ +£_)/2. For 
small displacements of the mechanical elements, Zf,x c <C 
£, the interaction Hamiltonian can be expanded to second 
order in x Cl yielding 



V 



Ho d f d c 
4tt £ 5 



-2Az 2 



12x 2 r + 180 



zjxl 
£ 2 



(3) 



where we ignored a constant term. The first two terms 
in Eq. §Ji§ describe frequency shifts due to the magnetic 
interaction. 

For high amplitude driving, the motion of the fork of 
effective mass mj and frequency ujf can be approximated 
by its classical amplitude. The cantilever motion, on the 



other hand, is treated quantum mechanically, with its 
displacement x c given in terms of the bosonic annihila- 
tion and creation operators b and 6^ by 



with 



x 
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where ui c and m c denote the cantilever's frequency and 
effective mass. The Hamiltonian governing the dynamics 
of the cantilever is then 



where 



H = Huj'M + V, 



lo c + Auj c 
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with Auj c the frequency shift from the magnetic interac- 
tion, 
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see Eq. In terms of b and w , the dipole interaction 
V becomes 



7T V 



where 



and 
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(10) 



(11) 



A final simplification is obtained by taking the driving 
frequency of the fork to be 



(12) 



then invoking the rotating wave approximation and 
switching to a frame rotating with w/, where V reduces 
to the familiar single-mode squeezing Hamiltonian (see 
e.g. M) 

V = h Xs {b^e~^ + b 2 e 12 ^ (13) 



with 
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and <p the relative phase between the classical fork and 
the cantilever. It defines which quadrature gets squeezed, 
and in the following we always consider the quadrature 
which experiences maximum squeezing, omitting the ex- 
plicit value of <f>. 
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B. Asymmetric case 

We now turn to the case where the equilibrium position 
of the cantilever nanomagnct is displaced by a distance 
£ from the center of the fork, the stability of that con- 
figuration being guaranteed as before by the mechanical 
stiffness of the cantilever. We now have 



£+ tzt + C + Zf, 

t- Xl£-(-Zf, 

and within the same limit as before 
Mo dfd c 
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When compared to the symmetric case, Eq. ([3]), the 
interaction (|1G[) comprises an additional squeezing con- 
tribution given by the term proportional to C z f x c- By 
driving at twice the cantilever frequency, and invoking 
the rotating wave approximation, one can access either 
of the two squeezing interactions separately, in a fashion 
reminiscent of the situation in parametrically coupled op- 
tomechanical resonators @. 

The resulting Hamiltonian is then the sum of two terms 
of the same form as in Eq. (|T3| with the coupling coeffi- 
cient of the term oscillating at 2uj' c given by 
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The coupling constant Xi(C) scales as the square of the 
amplitude of oscillations \f3\ 2 of the classical fork, rather 
than \j3\ as is the case for X2(£)j ancl hence can be domi- 
nant for strong fork driving. However, its dependence on 
£~ 7 rather than £~ 6 indicates that for appropriate values 
of £ the term proportional to \2 can be dominant instead, 
see Fig. 2. 



C. Experimental considerations 

Consider for concreteness a nanomechanical cantilever 
with natural frequency lu c = 2ir x 168 kHz, effective 
mass m c = 5 x 10~ 15 kg, and magnetic dipole moment 
d c = 2 x 1CP 14 A • m 2 . The fork is assumed to have ef- 
fective mass m/ = 1 x 10~ n kg, tips separation 2£ = 2.0 
/im, and magnetic dipole moments df = 1 x 10 -10 A • m 2 . 
These values result in a frequency shift of the cantilever 
of Aui c s=s 2tt x 72 MHz S> UJ C , indicating that the driv- 
ing frequency of the fork is dominated by the frequency 
shift ujf « Auj c . We further assume a fork oscillation 
amplitude of A f = 2z o \0\ = 10 nm, i.e. A f = O.OU, with 
associated mean phonon occupation |/3| 2 = 2.1 x 10 15 , 
which justifies a classical treatment. 




1 2 5 10 20 50 100 150 



FIG. 2. (Color Online) Coupling frequencies X^iO ( re(1 , 
solid), and Xi(C)(blue, dashed) in MHz as a function of dis- 
placement £ in nm for parameters of section II. C. 



For the asymmetric setup, we need to estimate the 
frequency of the cantilever in the ^-direction, which we 
obtain from u)c /oj c w w/h, where h and w are the 
cantilever's thickness and width. For h = 12 nm and 
w = 600nm this gives uj { c z) rts 2tt x 8.4MHz. With this 
value and Eqs. (fl"5j). we can find the critical points of the 
full interaction potential, Eq. (TTJ). For these parameters 
the maximum value of £ that yields a stable configura- 
tion is 162nm. The resulting strengths of the squeezing 
interaction are shown in Fig. [2] as a function of £• 



III. FLUCTUATIONS 

It is known that in the absence of fluctuations and 
starting from the oscillator ground state the squeezing 
Hamiltonian (|13l) produces a perfect "squeezed vacuum" 
with average phonon number 



<S + S> =sinh 2 (x*). 



(19) 



(In the following we use gcncrically the symbol x to de- 
scribe all setups of the previous section.) Clearly, this is 
unrealistic. What is missing from the discussion so far 
is a proper accounting of fluctuations. This section dis- 
cusses the effects of three sources of technical noise: the 
amplitude and phase fluctuations of the fork motion and 
the clamping noise resulting from the attachment of the 
cantilever to a thermal reservoir. 



A. Amplitude fluctuations 

We assume that the cantilever is slightly displaced from 
the center of the fork, driven at frequency ujf = 2uj' c , so 
that the squeezing strength is given by Eq. ([TT]) . The 
other two cases can be handled similarly and we give the 
results pertaining to them at the end of this subsection. 
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Assuming as usual that the amplitude of oscillations of 
the tuning fork can be decomposed into the sum of a 
constant amplitude and random amplitude fluctuations 



where 



\m) = \M+sm, 

we find the equation of motion 

i> = -[xo + Sx(t)} b\ 

where xo is given by Eq. (fl~7| and 

Sx(t) = nnW 
I Pol 
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Following Rcf. [Iff, we assume that the mean and two- 
time correlation functions of the amplitude fluctuations 
arc described by an Ornstein-Uehlenbeck process, so that 



(Sx(t)) = o, 
(Sx(t)Sx(t')) = ifViv 



r|i-t' 



(23) 
(24) 



where a is proportional to the variance of amplitude fluc- 
tuations and r is the inverse of the correlation time of 
the fluctuations. From Ref. [l6j . we find readily that the 
variances of the quadratures X\ and X 2 of the phonon 
mode are then given by 

(AA^i)) 2 = e- 2xat e- 2 K ^^(AA^O)) 2 , (25) 



(AA 2 (t)) 2 = e 2xo *e 2 /o d *'^(t')( A x 2 (0)) 2 



(26) 



which, upon carrying out the appropriate ensemble aver- 
ages and taking into account that the initial state of the 
cantilever is uncorrelated with the fluctuations, yields 
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Amplitude fluctuations reduce the rate at which the 
variances get squeezed. They do not however limit the 
maximum squeezing that can be generated. Moreover, 
Eq. (|29p shows that the lengthening in timescale result- 
ing from amplitude fluctuations scales as x§/|/3 | 2 , which 
is typically a small factor, indicating that in contrast to 
phase fluctuations to which we turn next, amplitude fluc- 
tuations in the classical drive do not significantly affect 
the dynamics of the system. Finally, we note that if the 
squeezing factor is given by Eq. (|18|) . we have to modify 
Eq. (|2"2"j) with an additional factor of 2, and the function 
f(t) of Eq. (j29j) is thus multiplied by a factor of 4. 

B. Phase fluctuations 



Again following Ref. [16] . we consider phase fluctua- 
tions 6<f>(t) about the relative phase <f>, approximated as 
a phase diffusion process characterized by the correlation 
functions 



(6</>(t)) = (<ty(i)> = 0, 
(6<f>(t)8<t>(t')) = D(t + t' -\t-t'\), 
(5<j>(t)5j>(t')) = 2D8(t-t'), 



(30) 



where D is the phase diffusion coefficient. From Ref. [17| . 
(see details in Appendix A) we find the variances of the 
quadratures 



(Ali(t)) 2 = - [cosh(2xi)exp(-L>t/2) - sinh (2 X t) exp(-5£>i/3)] , 
(AX 2 (t)) 2 = - [cosh (2 X t) exp(-L></2) + sinh (2 X t) exp(-5Dt/3)] , 



(31) 
(32) 



which reduce to 

(AAi(t)) 2 = ^Dtexp{2 X t) + ^cxp{-2 X t) 7 

(AA 2 (t)) 2 = icxp(2 x <) + ^Dtexp(-2 X t), (33) 

for Dt <C 1. Fig. depicts the squared quadrature 
(AA1) 2 as a function of the dimensionless time %t for 
several values of the diffusion coefficient D, illustrating 
the disappearance of squeezing for long enough times, an 
effect familiar from quantum optics. 



C. Clamping noise 

We now turn to a discussion of the effects of clamping 
noise on cantilever squeezing. We describe the thermal 
fluctuations coupled into the cantilever via a standard 
input-output formalism [l8j , resulting in the Heisenberg- 
Langevin equations 

S = -xSt-|S + Vo4>- (34) 
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FIG. 3. (AXi) 2 as a function of the dimensionless time yt for 
various amount of phase fluctuations: no fiuctuations(blue, 
solid); D = 10 _5 x ( green, dashed); D — 10 4 x (orange, dot- 
dashed); D = 10~ 3 x (red, dotted). 



Here 7 is the damping rate of the cantilever, and &i n 
a noise operator that accounts for thermal fluctuations, 
with 

(&*(*)> = $*(*)>=«, 

(SL(*)8ln(*')>=fith*(*-0. 
(Sin(*)6L(t')> = (*Hh + l)«(t-f). 
(^W^n(O) = <M*)MO> = 0, (35) 



and nth = kBT/fuv' c . This yields the quadrature evolu- 
tion equations 



d_ 

dt 1 
d_ 

dt 1 



^(AZO 2 = -(2x + 7)(A*i) 2 + 1 (2n th + 1) , (36) 



-^(AA 2 ) 2 = (2 X - 7)(AX 2 ) 2 + 7 (2nth + 1) , (37) 



and for a cantilever initially prepared in its ground state 
of ccnter-of-mass motion. 
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For t — > 00, (AXi) 2 approaches a steady state with 
reduced squeezing, as expected. Figure HI shows (AXi) 2 
for several thermal occupations of the clamp and two 
values of the mechanical damping rate 7. Consider for 
example a cantilever with a quality factor Q m = uj' c /2^ = 
10 4 , 7 ~ 2-7T x 3.6 kHz in a cryogenic environment at 
T = 20 mK. The phonon occupation number at ui' c is 
nth — 5.7. For a strong squeezing coupling constant of 



FIG. 4. (Color Online) (AXl) 2 as a function of scaled 
time with thermal fluctuations for different damping con- 
stants and occupation numbers: no fluctuations (blue, solid); 
fith = 5,7 = 10~ 2 x (green, dashed); n t h — 10, 7 = 3 x 10~ 2 x 
(orange, dot-dashed); and 7i t h = 10, 7 = 8 x 10 _2 x (red, dot- 
ted). 



X = 2ir x 2.8 MHz, this corresponds to a squeezing of 
2 x l(T 3 (21dB). 



IV. DETECTION 

The detection of motional squeezing could be per- 
formed along the same lines as in the experiments of 
Ref . , which generated classical squeezing in a paramet- 
rically driven mechanical cantilever and characterized it 
by measuring the two quadratures of oscillations using a 
fiber-optical interferometer. A more ambitious approach 
that offers many potential advantages involves a full 
determination of the cantilever state, rather than its co- 
variances only. One promising way to achieve this goal 
involves first transferring that state to an optical field, 
where detection techniques developed in quantum optics 
can then be applied. 

Quantum state transfer has already been the subject of 
a number of studies [l2l - fl4l . |20| . In particular, it is known 
that the two-mode beam-splitter interaction yields exact 
state transfer between two harmonic oscillators for ap- 
propriate interaction times and in the absence of dissipa- 
tion, and that this type of interaction can be realized in 
principle in the optomechanical interaction between the 
harmonically bound end-mirror of a Fabry-Perot and a 
near resonant intracavity light field. 

To implement that detection scheme, we expand the 
scheme of Fig. 1 to couple the cantilever to the intracav- 
ity light field of a resonator whose moving end-mirror is 
attached to the cantilever. This could be achieved e.g. 
with a cantilever of the type used in nanoscale magnetic 
resonance imaging see Fig. 5. Combined with a fixed 
large mirror, this can form a high finesse optomechanical 
resonator. Alternatively, one could also consider an ar- 
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FIG. 5. Schematics of the intracavity optical field optome- 
chanically coupled to the magnetically driven cantilever. 

rangcment where the cantilever forms a moving plate of a 
capacitor in a driven microelectronic LC-circuit [H, Q . In 
which case the coupling could actually be stronger, but 
one would be confronted with the lack of single-photon 
detectors in the microwave regime. 

In this section we consider two possible scenarios: In 
the first one the coupling to the optical field is present at 
all times, while in the second one a squeezed state of the 
cantilever is first prepared via magnetic dipole coupling 
to the classical tuning fork, and an optical field is sub- 
sequently turned on. We show that in both approaches 
the additional dissipation channel associated to the finite 
transmission of the Fabry-Perot leads to a significant re- 
duction in squeezing and imperfect state transfer. 



A. Continuous optical coupling 

With the additional optomechanical coupling of the 
cantilever to the optical field, the system Hamiltonian 
becomes 

H sys = H C + H p + H om + H m + H 1 +H K . (40) 

Here 

H c = huj a a^a (41) 

describes the optical cavity mode a of frequency loq, 

H p = ihr)(ah- iuJLt - ae iWL *) (42) 

accounts for the driving of the cavity by an external field 
of frequency usl at rate \q\ = ^Pn/Swi, with P the 
input power and k the cavity linewidth, and 

H om = -Kgtfa{b + P) (43) 



is the optomechanical coupling between the intracavity 
field and the cantilever mirror, with single-photon cou- 
pling frequency g. Finally, 7J 7 and H K describe the inter- 
action of the mirror and cavity field to thermal reservoirs 
and account for dissipation at rates 7 and k, respectively, 
and H m is the Hamiltonian of the magnetically driven 
cantilever of the previous sections. 

Considering a coherent pump of constant amplitude 
red-detuned from the cavity resonance by u' c , we expand 
the amplitudes of both the cavity field and cantilever 
oscillations as the sum of their expectation value and 
quantum fluctuations, 

a = (a) + Sd = Eq + 5a, 

b=(b) + So, (44) 

where (5a) = (5b) = and we neglect contributions from 
5a*5a compared to those from (a)5a) + (a)*Sa as usual, 
so that |-Eo| 2 = (o^){a) ~ (a^a). (Note that Eq is dimcn- 
sionless.) As is well known, this decomposition allows 
one to separate the optomechanical interaction Hamil- 
tonian into a classical Kerr-type contribution propor- 
tional to |i?o| 2 and a beam-splitter interaction, so that 
the cantilever dynamics is approximately described by 
the Hamiltonian 

H = -hg (EodaSP + E 6aUb\ + H K 

+ ^ (Sb^e-i* + S&e**} + H 7 , (45) 

where Eq accounts consistently for the Kerr nonlincarity. 
In steady state, the main consequence of the Kerr effect 
is a slight shift in the cavity resonance, an effect that 
can lead under appropriate conditions to optical bista- 
bility [2l|. Away from this multistable regime, the intra- 
cavity amplitude (a) = Eq of the radiation field is given 
approximately by (see Appendix B) 



k/2 — i[A + 2g 2 \Eo\ 2 /w' c ] ' v ' 

In the interaction picture, in the rotating wave approx- 
imation, and for uj' c — —A, the equations of motion for 
5a and 5b are, (see Eqs. (|B5[) and (|B6[0 

5a = — —6a + igEo5b + y/K&in, (47) 
5b = -l5b + igE*6a - X ^ + Vlhu, (48) 

where we have assumed that 

(al n (t)a in (t')) = 0, 

(M*)4,(*')> = *(*-* , )> 

(a in {t)a in (t')) = 0. (49) 

We consider for concreteness the case Eq = \Eq\c~ 171 ^ 2 , 
which holds when uj' c 3> K. The equations of motion for 
the variances of the position quadrature of the optical 
and phonon modes are then given by 

V = AV + B, (50) 
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with 
A 



B = 



V 



-k 2g\E \ 

-g\Eo\ _(k/2 + 7 /2+|x|) g\E \ 

-2g\E \ -(2|x|+7) 

k/4 


( 7 /8)(3n th + l) 



(AX 0>1 
(Al om , 

(AX m ,; 



(51) 
(52) 

(53) 



where (AX^i) 2 , i £ {o,om, m} are the variance of the 
position quadrature of the optical field, the covariance of 
position quadratures of the optical field and cantilever, 
and the variance of the position quadrature of the can- 
tilever, respectively. 

The steady-state squeezing of the intracavity field fol- 
lows from some elementary algebra, but its general form 
is cumbersome and we omit it here. In the physically 
relevant regime K 3> 7 it reduces to the simple form 



(AI M ) 2 (co) 



2r 2 



Is 1 



4(2r 2 + s)(l + 2s) : 



(54) 



where r = g\Eo\/K is the classically amplified optome- 
chanical coupling and s = x/ K - Figure [6] shows the 
steady-state squeezing of the optical field as a function 
of these parameters. It illustrates the monotonic increase 
in steady-state squeezing as s and r are increased, which 
is intuitively expected. However, this conclusion needs 
to be qualified by considering the steady-state limit of 
the cantilever squeezing. In the same limit it is likewise 
easily obtained as 



(AX m ,i) 2 (oo) 



(55) 



2(2r 2 + s)(l + 2s)' 

and is illustrated in Fig.[7]as a function of r and s. In con- 
trast to the situation for the optical field, we observe now 
that while increasing s increases the degree of squeezing, 
as would be expected, increasing the optomechanical cou- 
pling results in a decrease in squeezing. Indeed, while in 
the absence of optomechanical coupling the state of the 
cantilever mode would be almost perfectly squeezed, this 
ceases to be the case once the optical coupling is present. 
However, realistic experimental parameters yield enough 
squeezing transfer to the cavity field to be successfully 
detected. 

The physical origin of this behavior is that in addi- 
tion to the magnetic squeezing interaction the cantilever 
is now also subjected to the beam-splitter interaction. 
It results in a transfer of squeezing to the optical field, 
where it is now exposed to a dominating decoherence 
channel associated with the cavity loss rate n, normally 
much faster than the mechanical decay rate 7. This is 
more readily apparent in Fig. [8l which clearly illustrates 
how cavity damping decreases the steady-state squeezing 
of both the cavity field and cantilever for fixed beam- 
splitter and squeezing coupling constants g\E$\ and x- 
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FIG. 6. Steady-state squeezing of the cavity field as a function 
of the dimensionless coupling parameters r = g\Eo\/n and 
s = x/k. 
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FIG. 7. Steady state squeezing of the cantilever as a function 
of the dimensionless coupling parameters r = g\Eo\/n and 
s = x/k. 



B. Delayed detection 

The take-home message of the previous section is that 
while the coupling of the cantilever to the optical cavity 
allows detection of the squeezing, it does it at the cost of 
opening up a fast decoherence channel. For reasonable 
experimental parameters the resulting loss in squeezing 
is much larger than the limit imposed by thermal losses 
in the mechanics. This suggests that a better scenario 
might involve first preparing the cantilever in a strongly 
squeezed state, and only subsequently coupling it to the 
optical field. The issue with that approach is that it takes 
a time of the order of k _1 to switch on the intracavity 
optical field, a time during which the optical decoherence 
channel is already open. 

As before, we decompose the cantilever phonon field 
and intracavity optical field as the sum of their expecta- 
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FIG. 8. (Color Online) Steady state squeezing of the cav- 
ity field (red, solid) and cantilever (blue, dashed) as a func- 
tion of k/x hi the resolved-side band regime for fixed x an d 

g\E \/x = 9- 



tion value and quantum fluctuations, see Eqs. (|44[) . ex- 
cept that (a(t)) = Eo(t) is now an explicit function of 
time. The linearization process is questionable for very 
short times when the intracavity field is still extremely 
small. However, the optomechanical coupling is normally 
weak in that case, so that it should not qualitatively 
change the main features of the system dynamics. 

As shown in Appendix B, Eqs. (|B5|) and (|B6|) . the 
Heisenberg equations of motion for 8a(t) and 8b(t) are 
approximately given by 



8a = [iA - k/2] 8a + igE (t)(8b + Stf) + 



(56) 



6~b = [-ioj' c - 7 /2] 6b + ig [E£(t)5a + E {t)6a 1 ] 

- iixcos(u) f t + <j>){8b + 8P) + y/jb in , (57) 

and the evolution of i?o(i) is determined by Eqs. (|B12jl - 
(|B16|) . From these equations, it is possible to derive 
a closed set of equations for the first and second mo- 
ments of the operators 6a and 6b, see Appendix B. These 
equations could not be solved analytically, so this sub- 
section presents selected numerical results that illustrate 
the main features of the system dynamics. 

Fig. [9j which is for a relatively high-loss optical cavity 
that allows for a fast switching of the optical field, shows 
the coupled dynamics of the cantilever and optical fields 
in a situation where the cantilever was first prepared in 
a squeezed state, before the optical field is switched on 
at = X ■ It illustrates a situation where squeezing 
transfer suffers from the broad decoherence channel of the 
optical cavity. Thus the squeezing is not efficiently de- 
tectable in the cavity field. Note also that since the beam- 
splitter interaction frequency g\Eo\ — n in that example, 
the oscillatory coherent state transfer between the can- 
tilever and the optical field is strongly suppressed, with 
the energy of the cantilever-field system being rapidly 
lost through the optical decay channel. A much more 
significant coherent exchange between the two subsys- 
tems requires cither a stronger field amplitude |£b|i or a 




FIG. 9. (Color Online) Squared quadratures of position for 
the cavity field (red, solid) and the mechanical oscillator (blue, 
dashed) as a function of scaled time. The optomechanical 
coupling is turned on at the dimensionless time = 1 and 
coherently builds up towards g\Eo\/x = 9, k/x = 10. 



slower decay of the light field, so that g\E \ ^> k. Such an 
example is illustrated in Fig. 1101 which shows the charac- 
teristic coherent state transfer between the phonon and 
photon field, as expected. One problem here is of course 
that by decreasing the cavity damping rate, one requires 
a longer time to turn on the light field to its final value 
\Eo\, thereby increasing the role of dissipation. Still, in 
this situation it is possible to achieve a reasonably good 
transfer of squeezing from the cantilever to the optical 
field. 

During the coherent state transfer, the maximum 
squeezing in the intracavity field occurs after half an ex- 
change period. In Fig. 111! the minimum values of the 
quadrature variance in the cavity field are plotted as a 
function of g\Eo\/n. As is expected, smaller cavity damp- 
ing rate or lager coupling strength gives rise to stronger 
maximum squeezing in the cavity field. 



V. CONCLUSION 

In summary, we have presented a theoretical analysis 
of the motional squeezing of a cantilever magnetically 
coupled to a classical tuning fork via microscopic mag- 
netic dipoles. We showed that this coupling can result 
in significant squeezing of a quadrature of motion of the 
cantilever if appropriately driven by a classical force, and 
found that the system is robust against various sources of 
noise, with phase noise in the driving of the classical driv- 
ing tuning fork the dominant source of decoherence. We 
proposed a scheme for the detection of the effect based 
on state transfer to the intracavity field of an optical 
resonator with one end-mirror formed by the oscillating 
cantilever. Challenges to the measurement process asso- 
ciated with the additional decoherence channel opened 
by the coupling to the optical resonator were discussed. 
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FIG. 10. (Color Online) Squared quadratures of position in 
the cavity field (red, solid) and the mechanical oscillator (blue 
dashed) with strong coherent optomechanical coupling. Here 
<?|£ |/X = 9,«/X = 1- 



chanical coupling to control and probe the system, in- 
cluding the use of multimode light fields. 
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Appendix A: Effect of phase fluctuations 

This appendix presents details of the evaluation of the 
effect of phase fluctuations for the case of an asymmetric 
setup. The symmetric case is analogous. In the pres- 
ence of phase fluctuations, the Hcisenberg evolution of 
the phonon annihilation operator b(t) becomes 




g\E \/K 



FIG. 11. Minimum values of squared quadrature variance in 
the cavity field plotted as a function of g\Eo\/K. 



It has recently been proposed that pulsed optomechan- 
ical configurations permit mapping the quantum state of 
optomechanical oscillators by using a sequence of appro- 
priately shaped optical pulses separated in time by half 
a vibration period of the mechanical system (22|. This 
approach presents several advantages, the first one being 
that it allows for the use of low finesse optical resonators 
that permit fast switching and the second being that the 
oscillator is coupled to the optical dissipation channel for 
very short times only. Unfortunately, this scheme relies 
on the mechanical oscillator being subject only to free 
evolution between the light pulses, which is not the case 
here since the squeezing interaction is acting at all times. 
It is not clear how it could be rapidly be switched off 
to make the pulsed detection scheme applicable. Future 
studies will consider whether adapting pulsed detection 
scenarios to the present situation may be possible. We 
will also consider the quantum dynamics of the tuning 
fork, the interaction between these two subsystems on a 
quantum level as well as further possibilities of optome- 



dt 



b(t) 



iS<t>{t] 



ftt( t ). 



(Al) 



and the full dynamics of the system can be expressed in 
the general form 



f = (A + i5j>(t)B)Y, 



(A2) 



where Y is the vector of bilinear operators I r, w 2 , and 
b'b + bb* . Due to the form of the Heisenberg equations 
of motion for the problem at hand, the phase fluctua- 
tions 5<j){t) can only be eliminated in one bilinear entry 
in Y, and it is therefore impossible to readily perform 
the statistical average over phase noise. As shown by 
Wodkievicz [ijj this difficulty can be circumvented by 
solving two systems of matrix equations separately for 
two particular choices of Y, b and wb + bw . For b 2 we 
have 



Y = 



A 



e -2i5<t>{t)y\2 



-x o 
-2 X -2 X 
-x o 



B 




-1 




(A3) 



(A4) 



(A5) 



while for b^b + btf 



Y 



6&t + b% 

e is<t>{t)yi 

e -i50(t)Jf2 



(A6) 
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.4 



B 



-2 X -2.x 
-X 
-X 




1 
0-1 



(A7) 



(A8) 



If the initial state is uncorrelated with the fluctuations, 
as is physically the case, this form of equations can be 
solved exactly to give 



(Y(t)) 



,(A-DB 2 )t 



<no)>, 



so that 



(6 2 ) ps --sinh(2xt)exp[-5Dt/3] 



(A9) 



(A10) 



(S f 6 + bti) PS cosh (2 X t) exp[-Dt/2], (All) 
where we have assumed that D -C X f° r simplicity. 

Appendix B: Equations of motion for optically 
coupled system 

In the rotating frame at frequency uj, the system 
Hamiltonian (|40p is 

H = -hAa)a - hga)a(b + P) + ifvq(a) - a) + flu^S 
+ /iXCOs(w / i + 0)(S + S t ) 2 + H 7 + i? K , (Bl) 
where 

A = uj l - oj , (B2) 
and the equations of motion for a and b are 

a = iAa + ig(6 + P)a — —a + 7] + s/Kain, (B3) 



b = —iio' c b + iga'a 



- 2i X cos(oj f t + (f))(b + P) - ^b + ^fb in , (B4) 

where <£> is the phase difference between the tuning fork 
and cantilever. Introducing 5a = a — (a) and Sb = b — 
(b) yields then for 5a and 5b the linearized equations of 
motion 

5a ss [iA - k/2] 5a + ig(b + P)5a 

+ ig(a) (5b + 5P) + \fKa- m (B5) 
Ps [iA - k/2] 5a + ig(a) (Sb + 5P) + v^ain, 

Sb m [-iw^ - 7/2] 56 + ig((a)*5a + (a)5a?) 

- 2i X cos(ujft + (t))(5b + 5P) + ^yb in . (B6) 



The second, approximate form of Eq. (|B5[) results from 
the fact that the mean phonon number is of order unity, 



which is much smaller than the mean number of intrac- 
tivity photons (10 4 — 10 8 ). Under these conditions one 
can neglect the term ig(b + 6 T ) in that equation. 

From these approximate linearized equations we easily 
obtain the equations of motion for expectation values of 
the quadrature operators of a and 6, 



(X a ) = -A(Y a ) - -(X a ) - 2g(X b Y a ) + r h (B7) 
(Y a ) = A(X a ) - ^(Y a ) + 2g(X b X a ), (B8) 

{£ b ) = w' c <y b )-l{x b ), (B9) 

(A) =~K + 4xcos( W/ i + 0)] (X b ) 

+ g(N a )-l(Y b ), (B10) 

(N a ) = 2 V (X a ) - K {N a ), (Bll) 

where X a = (a + a')/2, Y a = (a — a')/2i, and N a is 
the intracavity photon number operator. In the regime 
where uj' c » x S> 7 and (N a ) = \E Q \ 2 > 1, the change 
in the classical component of the intracavity field due 
the cantilever oscillations remains very small. It can be 
ignored in determining the dynamics of the beam-splitter 
coupling constant g\Eo(t)\, which is then governed by the 
approximate equations of motion 



(X a ) = -A(Y a ) - -(X a ) - 2g(X b Y a ) + v , (B12) 
(Y a ) = A(X a ) - |(F a ) + 2g(X b X a ), (B13) 



(X b )=oj' c (Y b ), 
(Y b ) = -Lu' c (X b )+g\E \ 2 , 
\E \ 2 = 2ri(X a )-K\Eo\ 2 . 



(B14) 
(B15) 
(B16) 



These equations yield the time evolution of (a(t)), as well 
as its classical steady-state value. We find, upon factor- 
izing (X b Y a ), 



Eq\ss — 



n 



n/2-i[A cfi + 2g 2 (\E \l)/u;' c ] 



(B17) 



In the rotating frame at frequency lul, the position 
quadratures of both cantilever and cavity field arc defined 



-iAt 



Ai j0 = -^(dae 
X hm = 1 (5be 1 ^ 



5ate iAt ), 



SPe'^ 1 



), 



(B18) 
(B19) 



In order to calculate their variances we need to have the 
expectation value of second moments of the fluctuations. 
Taking quantum averages of the equations of motion for 
these quantities results in the closed set of equations 
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— (M 2 ) = 2iA(Sa 2 ) - k(8ci 2 ) + 2ig(a)(SaSb + SaStf), 

— (SalSa) = -k(5o) So) - ig \ (a)* (SaSb + SaSP) - (a)(Sa i Sb + SaUP) 
at ' L 



dt 

d 
dt 



-2iw' c (5b 2 ) - 1 (5b 2 ) + 2ig Ua)*(5aSb) + (a)(Ja t <SS) 
2i x cos(w/t + <j)){25b 2 + SPSb + SbSP), 



— (SPSb) = -~f(SPSb) - ig (a)* (SaSb - SaSP) + (a) (5tf 5b - Stfstfl) 
dt L 

+ 2ixcos(uj f t + (j))(Sb 2 - Sb t2 ) + -fn th , 
A (SaSb) = i(A-cj' c ){5aSb) - [(k + 7 )/2] (SaSb) + ig \(a)(5b 2 + 5PSb) + (a)*(Sa 2 ) + (d)(8a8a^ 

—2ix cos(w/t + 4>){5d5b + SaSb'}, 
(SaStt) - t(A + u' c ){8a8fr) - [(« + 7 )/2] (SdSP) + ig \(a)(SP 2 + SW) - (a)* {6a 2 ) - (a)(5aSa^) 
+ 2i\ cos(w ft + 4>) (SaSb + SdSP) , 



(B20) 
(B21) 

(B22) 
(B23) 
(B24) 



(B25) 
(B26) 



r 



These are the equations that are solved numerically to obtain the figures of section IV. B. 
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